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Abstract 

Self-consistent theory of electron localization in disordered systems is gen- 
eralized for the case of interacting electrons. We propose and critically com- 
pare a number of possible self-consistency schemes which take into account the 
lowest perturbation theory contributions over the interaction. Depending on 
self-consistency scheme we can obtain either the continuous metal-insulator 
transition or that with the minimal metallic conductivity. Within the continu- 
ous transition approach we calculate the frequency dependence of generalized 
diffusion coefficient both for metallic and insulating phases. We also con- 
sider interaction renormalization of the single-electron density of states which 
demonstrates the growth of the effective pseudogap as system goes from metal 
to insulator. 

PACS numbers: 71.30.+h, 71.55.Jv, 72.15.Rn 
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I. INTRODUCTION 



Theoretical description of disorder-induced metal-insulator transitions constitutes one of 
the major problems of condensed matter theory [1,2]. The main difficulties here are con- 
nected with consistent account of interelectron interactions, the major importance of which 
was already demonstrated even in case of weakly disordered metals [3] . In recent years this 
problem was actively studied within the renormalization group approach [1,2,4], generalizing 
the usual scaling theory of localization [1]. Though very successful in many respects this 
approach was not able to produce a complete solution of the problem and in particular it 
almost failed in its ability to predict the concrete and experimentally verifiable dependences 
of physical properties on parameters of the system which control the transition. Renormal- 
ization group (scaling) approach is well fitted for the analysis of the region very close to 
transition itself [2], while experimentally we are usually dealing with variation of physical 
properties in some wide region of parameters controlling the transition. In particular, up to 
now there are no papers treating the insulating phase within renormalization group approach 
accounting for interactions. 

In the theory of metal- insulator transitions neglecting the interactions along with scaling 
approach [1], the so called self-consistent theory of locahzation is widely used [5,6,7] and pro- 
vides a rather successful interpolation scheme, which allows calculations of the main physical 
properties of a system for the wide region of parameters from weakly disordered metal to the 
Anderson insulator. This theory also reproduces all the major results of the scaling approach 
[5,6]. While the rigorous diagrammatic procedure leading to the self-consistent scheme is 
still unknown, this approach leads to results which are in quantitative accordance with ex- 
act numerical data modelling the Anderson transition [8,9]. First attempts to account for 
the interaction effects within the self-consistent theory of localization were undertaken in 
Refs.[10,ll] (Cf. also Ref [6]). 

This paper attempts to formulate some version of self-consistent theory of disorder- 
induced metal-insulator transition accounting for the first-order perturbation theory correc- 
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tions over interelectron interactions. The difference with the previous attempts [10,11] is that 
we shall try to take into account the interaction effects upon the generalized diffusion coef- 
ficient, which represents the main physical characteristic to be determined self-consistently 
Unfortunately, we can propose a number of alternative formulations of self-consistency which 
lead to different predictions concerning the physical properties. We can only choose between 
these different schemes only using some guidance from the experimental behavior during 
metal-insulator transition. In particular, depending upon the self-consistency scheme we 
can get either the continuous metal-insulator transition or that with the so called minimal 
metallic conductivity. Modern experiments lead to the picture of continuous transition [1,2]. 
In the framework of such a scheme we shall calculate the frequency dependence of generalized 
diffusion coefficient as well as renormalization of single-electron (tunneling) density of states 
due to interactions for the wide region of parameters varying during the metal-insulator 
transition. 

II. INTERACTION CORRECTIONS TO CONDUCTIVITY CLOSE TO THE 

ANDERSON TRANSITION. 

Let us introduce two-particle Green's function for electrons in the random field of impu- 
rities [5]: 



Let us introduce the appropriate full vertex-part rpp,(qa;) and also the "triangular" 
vertex (pqo;) (Cf. Fig.l). 

Within self-consistent theory of localization the full vertex-part and triangular vertex 
are given by [6]: 




(1) 



where p± = p ± |, and <...>- denotes the averaging over the random potential. 




(2) 
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where 1 — — TrpU^No{0) - is the usual Born decay rate, p - impurity concentration, U - 
impurity potential, No{0) - density of states at the Fermi level for noninteracting electrons, 
D{u!) - the generalized diffusion coefficient. 

Expressions for the vertex parts actually coincide with those obtained in the "ladder" 
approximation with the usual Drude-like diffusion coefficient Dq replaced by frequency- 
dependent generalized diffusion coefficient D{uj), defined from the solution of the following 
self-consistent equation [5,7]: 

D{u) 7r«-„(0) ^^^^ -'m + D{,^)q' ^ ' 

where — 'min{l~^, kp} - is integration cutoff momentum, I - is the mean-free path, kp - 
is Fermi momentum. 

In three-dimensional system the generalized diffusion coefficient takes the following form: 

^(f^Y^a uj<^ujc, O(>0 Metal 

LJ ^ LJc Metal and Insulator (5) 

^(l^)"' (^) ^ < ^c, Oi<0 Insulator 

where A = - is dimensionless disorder parameter, Ep - Fermi energy, cj^ = " 
the characteristic frequency, a — 1 — SXxq - is the parameter which controls metal-insulator 
transition, xq — ^ - the dimensionless cutoff. 

Neglecting localization contributions first-order corrections to conductivity due to 
Coulomb interaction are determined by diagrams shown in Fig. 2 [12]. It was shown in 
Ref.[12] that the total contribution of diagrams (a), (b) and (c) is actually zero and con- 
ductivity correction reduces to that determined by diagrams (d) and (e). Here we neglect 
also the so called Hartree corrections to conductivity [3,12], which is valid in the limit of 
2kF/K ^ 1, where k - is the inverse screening length. This inequality, strictly speaking, is 
valid for systems with low electronic density, which are most interesting for experimental 
studies of disorder induced metal-insulator transitions. Also, if we remember the known 



Do 



2 \Ef J 
Ef 



2 



results on the divergence of screening length at the metal-insulator transition, we can guess 
that this approximation becomes better as we approach the transition. The point-like in- 
teraction model used below has to be understood in this sense. 

Using the explicit form of impurity vertexes and we obtain the following correction 
to conductivity due to interactions (Cf. [12]): 

^""^"^^ - I 2^1 {27^Y {-liSl + uj) + DiSl + u)q^){-in + D(^])g2)2 

Here e - is the electronic charge, c? - is spatial dimensionality, V{qVL) - interaction. For 
simplicity in the following we shall mainly consider the point-like interaction V{qVL) = Vq. 
In case of dynamically screened Coulomb interaction at small q and Q we have [3]: 

Then for the first-order correction to diffusion from the point-like interactions we obtain: 

5Dc{uj) 5a{ijj) 



Do 2emoiO)Do 

1 f rr^ f d'^Q <f 



(8) 



oo 



Tid'^ ttNo{0)J J {2ttY + u) + D{n + uj)q^){-in + D{n)q^y 

UJ 

where /i = VoA^o(O). 

Using in (§) the form of diffusion coefficient obtained in self-consistent theory of local- 
ization flSD in three-dimensional case we obtain: 



Do ^3ttX 



3V^±2V2fE^Y c,>0 Metal 

^(-^Y uj:^uJc Metal and Insulator (9) 

^<^c, «<0 Insulator 

It is easy to see that the correction to diffusion coefficient (conductivity) diverges as system 
approaches the Anderson transition from metallic phase —>■ 0), while in insulating phase 
divergence appears as — > 0. 

If instead of the point-like interaction we take dynamically screened Coulomb interaction 
(0), correction to diffusion coefficient takes the form: 



oo feo 



Do ^ TT kp J ^^I ^'^{-i{n + to) + D{n + uj)q^){-tn + D{n)q'')D{n) 

uj 

Using here the form of diffusion coefficient from self-consistent theory locahzation (|^), we get 
the correction hke (|]), but with /i = |. That is, the use of dynamically screened Coulomb 
interaction instead of point-like leads just to a replacement of /i by a constant of order of 
unity and in the following we always assume the point-like interaction. 

Thus we can see that the correction to diffusion coefficient close to the Anderson transi- 
tion becomes much larger than diffusion coefficient of the self-consistent theory itself. This 
obviously leads to the necessity to formulate some kind of new self-consistency scheme taking 
interelectron interaction into account from the very beginning. 

III. SELF-CONSISTENCY SCHEMES. 

Let us start with the usual self-consistent theory of localization in the absence of Coulomb 
interaction [5 — 7] . In case of weak disorder conductivity is determined by Drude diffusion co- 
efficient Dq. Summation of "maximally crossed" diagrams leads to the following localization 
correction to diffusion coefficient [5]: 

Do 7rNo{0)^^,^^-tuj + Doq^ ^ ^ 

Now we can introduce the so called relaxation kernel M{uj), which is connected with the 
generalized diffusion coefficient by the following relation [5]: 

In particular, the Drude relaxation kernel is given by Mq = = 2^7. 

Correction to relaxation kernel can be expressed via the correction to diffusion coefficient 



as: 
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Consider the usual Drude metal as the zeroth approximation: 



5M{uj) = -^6D{uj) (14) 



Replacing Drude diffusion coefficient Dq in the diffusion pole of Eq. (|TT]) by the generalized 
one D{uj), and using this relation in Eq.(|T^, we obtain the main equation of self-consistent 
theory of localization: 

which coincides with Eq.(|) if we take into account that ^^^^ = 77^- 

The first order corrections due to Coulomb interaction are determined by diagrams 
(d) and (e) shown in Fig.2. Unfortunately we cannot propose the uniquely defined self- 
consistency scheme which take into account corrections due to interelectron interactions. 
We can only formulate several alternative schemes and try to choose between them using 
some additional qualitative guesses. 

Scheme I. Let us first use the approach similar to the usual self-consistent theory of 
localization. As a zeroth approximation we take Drude metal and consider localization and 
Coulomb corrections on equal footing and introducing the self-consistency condition in terms 
of relaxation kernel. 

In this case the relaxation kernel takes the following form: 

M{uj) = Mo + 6M{uj) (16) 

where SM{uj) = SMi{uj) + SMc{uj) = -^{SDi{uj) + SDc{uj)). Here the localization correction 
to diffusion coefficient Diioj) is defined by Eq.(|l^), while the Coulomb correction Ddyj) is 
given by Eq.(P). Self-consistency procedure is reduced to the replacement of Dq by the 
generalized diffusion coefficient in all diffusion denominators. As a result we obtain the 
following integral equation for the generalized diffusion equation: 

-5^ = 1 + ^-/^ \ (17) 

D{^) iriV|,(0) / {2ttY -iu] + D{^)q' ^ ' 
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TTd"^ nNo{0)J J {2TTY{-i{n + uj) + D{n + uj)q^){-tn + D{n)q^y 
Scheme II. As a zeroth approximation we can now take the "dirty" metal described by 
the usual self-consistent theory of localization and consider first-order Coulomb corrections 
to this state. Self-consistency again is realized through corrections to relaxation kernel. 

Usual self-consistent theory is based upon Eq.([TB|). Now we have to add into the right- 
hand side of this equation the Coulomb correction to relaxation kernel 5Mc{uj). This cor- 
rection has the following form: 

5M,{u:) = -^f^SD^ico) (18) 
L>s[UJ) 

where Ms{uj) and Ds{uj) - are the relaxation kernel and diffusion coefficient of the usual 
self-consistent theory of localization (^. 

In this case, the equation determining for diffusion coefficient reduces to: 

Do 1 f d'^q 1 



D{uj) 7riVo(0) J {27fY -luj + D{uj)q 

Do Y 8i ^ 1 1 r d'^q q 



\Ds{uj) Tid"^ ^7tNo{0)J J {27rY + uj) + D{Q + uj)q'^){-iQ + D{Q)q'^y 

Scheme III. Now we can treat the Coulomb correction to relaxation kernel 6MJuj) in 



the right-hand side of self-consistency equation ([ISD also in self-consistent manner 



SM^iu) = -MMsD^iu;) (20) 
U[UJ) 

In this case the equation for diffusion coefficient becomes: 

^° 1 + / fracd'q{2nr - (21) 



D{uj) ttNo{0)J ' ' -iiu + D{iu)q^ 

Do V St ^ 1 f d^q g2 



— -mDq — / d^l 



\D{uj) TTd"^ '^ttNo{0)J J {2TTY{-i{n + uj) + D{n + uj)q^){-in + D{n)qY 
Scheme IV. As a zeroth approximation we can again take the usual self-consistent 
theory of localization but self-consistency scheme can be realized considering corrections to 
diffusion coefficient itself, not to the relaxation kernel. 
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Self-consistent equation (H) can be rewritten as: 

D(uj) 1 



(22) 



Da l-J ^— f 1 

" ^ 7rAfo{0) J {2ttY -iuj+D(u))q^ 

Let us now add to the right-hand side the correction to diffusion coefficient Eq.(§), so that 
finally we obtain the following equation: 

1 + (23) 



^0 1 + 



tA^o(O) ■> {2-kY -iuJ+D{uj)q 



oo 



The choice between alternative schemes I-IV is difficult to make on general grounds. 
Thus, first of all, we have to analyze these schemes qualitatively. 

IV. QUALITATIVE ANALYSIS OF DIFFERENT SELF-CONSISTENCY 
SCHEMES IN METALLIC PHASE. 

In the following we consider only three-dimensional systems. In self-consistent theory of 
localization the diffusion coefficient is determined by Eq.^). From here it is easy to find the 
relaxation kernel: 



Ms{uj) Do 



(24) 



Mo Ds[uj) 

Consider now the correction to relaxation kernel which is given by the following expression: 

^Mo(cu) ^ 5D,{uj) 

Mo Do ^ ^ 

(As a zeroth approximation we take the usual Drude metal here.) 

The Coulomb correction to diffusion coefficient of the usual self-consistent theory of 
localization is determined by Eq.(^. Using this in Eq. (^), we obtain the appropriate 
correction to relaxation kernel 6Mos{uj). 

Comparing 6Mos{uj) with Ms{uj), it is easy to see that in metallic region and for the 
whole interval of frequencies we have the following qualitative estimate: 



SMosiuj) ^ cfiMsiuj) (26) 

where c - is a numerical coefficient of the order of unity. 

Detailed analysis shows that the similar estimate in metallic phase is valid for relaxation 
kernel correction in all schemes discussed above. The reason here is that all schemes of 
self-consistency in metallic phase the frequency behavior of diffusion coefficient is similar to 
that of the usual self-consistent theory of localization — it is a constant for small frequencies 
and D{uj) ~ (— iu;)^/^ for high-frequency region. Thus, instead of Eq. (p6D we can write: 

SMo{uj) ^ ciiM{uj) (27) 

Using the expression (pTf ) we can analyze the qualitative behavior of different self-consistency 
schemes in metallic phase for = 0. 

Scheme I. We are considering the metallic phase, thus for = we have D{uj = Q) = 
D^O. Muhiplying Eq.{^) by Mq, we obtain: 

Mo f d'^q 1 



M = Mo + — ^ /-^A + <5Mo(0) 



7rA/'o(0) J {2TiYDq^ 
Using Eq.(l^), we get: 

M = Mo + SAxoM + c/iM 

Accordingly M = l^^xxo-c^l 

— = 1 — 3Aa;o — c/x = a — a* (28) 

where a* = c/x. 

Thus, even in the presence of weak Coulomb interaction continuous Anderson transition 
persists and conductivity exponent remains the usual v = 1 (i.e. diffusion coefficient linearly 
goes to zero with disorder parameter a — a*). However, we observe the shift of transition 
into the region of more weak disorder a = a* . Dependence of diffusion coefficient on disorder 
is shown in Fig. 3 (a). 

Note that if we analyze Eq. (p!7| ) assuming the generalized diffusion coefficient independent 
of frequency, we shall obtain discontinuous metal-insulator transition with minimal metallic 
conductivity. 

10 



Scheme II. Multiplying Eq.(|^) by Mq and using Eq.(]27|) we obtain: 

M = Mo + SAsqM + ^c/iM 



From here we have M = — or 



D cjj, 



n -^--2 (29) 



In this case again the continuous Anderson transition persists, conductivity exponent is 
again u = 1 and we only have the shift of the transition to the region of weaker disorder 
a = (c/i)i/3. 

Dependence of diffusion coefficient on disorder is shown in Fig. 3(b). 
Scheme III. Multiplying Eg. (pi]) by Mq and using Eq.(^), we obtain: 

M = Mo + 3Aa;oM + hrp c/xM 

or 

M /M\^ 
Mo ^\mJ 

In terms of diffusion coefficient: 

a^— + (30) 



In this case we get the minimal metallic conductivity: 

^ = (2c/i)i/=' for a = «* = ^(2c/i)^/3 

Dependence of diffusion coefficient on disorder is shown in Fig. 3(c). 

Scheme IV. Diffusion coefficient correction is obtained from Eq.(p5D using Eg. (pT]) : 

SDcjuj) ^ 6Mo{lu) ^ cfiMju) ^ _ _Do_ 
Do Mo Mo ^^D{u) ^ ^ 

Substituting (|3) into Eq.(|2|) we get: 



D ^ 1 Do 

Do l + {l-a)^ ""^D 
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or 



" = TTTf + 7^ (32) 



C/i 

Again we obtain the minimal metallic conductivity: 

^ ^ (2c/i)i/3 for a = «* ^ ^(20/1)^/=^ 
Dependence of diffusion coefficient on disorder is shown in Fig. 3(d). 

Qualitatively the results of schemes III and IV coincide and lead to the concept of 
minimal metallic conductivity which does not correspond to the majority of experimental 
data we have in this field [1, 2]. Schemes I and II lead to the same continuous behavior of 
diffusion coefficient, differences appear only in corresponding estimates of critical disorder 
at the transition. It is quite natural that our approximations lead to the increased tendency 
toward metal-insulator transition due to interelectron interactions. The transition takes 
place at smaller disorders than in the absence of interaction. 

In the following we shall concentrate on self-consistency scheme I where Coulomb and 
localization corrections are being treated on the equal footing. 

V. FREQUENCY DEPENDENCE OF DIFFUSION COEFFICIENT. 

In self-consistency scheme I diffusion coefficient is defined by the integral equation (p!7D. 
Let us transform it to dimensionless Matsubara frequencies: cu, Vl, and also 

introduce the dimensionless diffusion coefficient d{uj) = In these notations integral 

equation ([T7|) takes the following form: 

1 } dyy'-' , 



1 + —-dXx 



8 x d-2 7 / dyy'^^' . 

+ ° J d{u + Q)d^{Q)J f,,2 , A / 2 , n ^ ^ ^ 
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In the following we shall limit ourselves only to the case of spatial dimension d=3. Diffusion 
coefficient of the usual self-consistent theory of localization (|]) in the same notations reduces 
to: 

a = 1 — 3Axo oj <^ ojc, a > Metal 

2 

(^|3Aa;o) ^ 3> cUc Metal and Insulator (34) 

(■^3Aa:o) ^ _ j^^^^^2^ UJ ^ UJc, « < lusulator 



d{uj) 



where Uc = ,J"^ s-i and ^ - is the localization length. Let us introduce K(uj) = -rf-r and 



, 2 



analyze Eq.(^) assuming that K{uj), K{Q) and K{uj + Q) <^ 1. Expanding the right-hand 
side of Eq. (Bl) over these small parameters we obtain: 



d{u) 2 d{Lu) 

+ ^^^"^7 d{u + Q)d^{Q) (irV2(n) + Ry^iuj + n)f ^ 
Consider the metallic phase and look for diffusion coefficient d{u!) solution in the following 
form: 

{d U <^ LUc 

/ ,1 (36) 

Substituting (^) into Eq. (|35|) we find d and cUc and for the diffusion coefficient we obtain: 

a — a* uj <^ CO.. 

d{u;) = { , 2 ^ (37) 



2 

where Uc = U , a* = cfi, c ~ 0, 89 

(f 3Azo) 

Thus for the metallic phase we confirmed the previous qualitative conclusion — Anderson 
transition persists and the conductivity exponent remains u = 1. The transition itself has 
shifted to the region of weaker disorder a = a* = cfi. The frequency behavior of diffusion 
coefficient in metallic phase is qualitatively similar to that in the usual self-consistent theory 
of localization (P^. In the region of high frequencies u ^ Uc the behavior of diffusion 
coefficient remains unchanged after the introduction of interelectron interactions. 
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Consider now the insulating phase. In the region of high frequencies uj ^ uOc the diffusion 
coefficient obviously possess the frequency dependence like d{(jj) ~ uj'^^^. Assume that for 
small frequencies it is also some power of the frequency: 

(38) 

where S is some exponent to be determined. 

Substituting ( p8D into (|35D and considering the case of a < (insulating phase of the 



usual self-consistent theory of localization) and |a| a*, we get: 



d{u;) = { - . (39) 

I 1 3 2 

where Uc = , .i , while uj* ^ 0, lu-. — - — = 0, Ittt^ — is some new characteristic 
frequency defined by the interactions. Note that cj* — * as we approach the transition 
point when ^ ^ oo. 

Thus, sufficiently deep inside the insulating phase when a < |a| ^ a* and for the 
frequencies u ^ u*, the diffusion coefficient remains the same as in the self-consistent 
theory of localization, i.e. at small frequencies it is linear over frequency, while for the 
higher frequencies it is ~ u^^^. 

The analysis of Eq.(p5|) shows that for the frequencies uj <^ u* it is impossible to find 



the power-like dependence for d{uj), i.e. the diffusion coefficient in the insulating phase is 
apparently can not be represented in the form of d{uj) = d^(^-^^ , where 5 - is some 
unknown exponent. Because of this we were unable to find any analytical treatment of 
Eq.(p5D in the region of u u* within the insulating phase. 

Consider now the system behavior not very deep inside the insulating phase when a — 
a* < while a > 0, that is when the system without interaction would be within the 
metallic phase. 

Let us assume that the frequency behavior of the diffusion coefficient for u Uc possess 



the power-like form, i.e. the diffusion coefficient is defined by the expression (|38|) . 
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Substituting (|38D into (^) we get <5 = |- As a result for the diffusion coefficient we get: 

d{u) = { ^/ \l , (40) 



where uOc = -7 r^. 

(fSAxo) 

Naturally, the exact solution for the diffusion coefficient should show the continuous 
change of frequency behavior around uo ~ uOf.. 

Thus, within the insulating phase close to transition point, where the system without 
interactions should have been metallic, the diffusion coefficient behaves as ~ u^/^^ every- 
where, though for the low frequency region the coefficient of u;^/^ differs from that of the usual 
self-consistent theory of localization and explicitly depends upon the interaction constant. 

Note that if we use the dynamically screened Coulomb interaction when ~ 1, the 
region of applicability of (|40|) widens, because of a* ~ /i ~ 1. 

We have also performed the numerical analysis of the integral equation (|33D for the wide 
region of frequencies, both for metallic (Fig. 4) and insulating phases (fig.5). Numerical data 
are in good correspondence with our analytical estimates. 

In the region of high frequencies, both for metallic and insulating phases, the frequency 
behavior of diffusion coefficient is very close to that defined by the usual self-consistent 
theory of localization. 

In the region of small frequencies within the metallic phase diffusion coefficient d{uj) 
diminishes as interaction grows. Dependence of static generalized diffusion coefficient on 
disorder for /i = 0, 24 is shown at the insert of Fig.4, and is practically linear. Metal- 
insulator transition in this case is observed at a = a* = /i, where c ~ 0, 5, which is also in 
good correspondence with our qualitative analysis. 

Within the insulating phase for the region of small frequencies {uj <^ uj*) we observe 
significant deviations from predictions of the usual self-consistent theory of localization. 
Diffusion coefficient is apparently nonanalytic in frequency here and we clearly see the ten- 
dency to formation of some kind of effective gap for the frequencies uj <^ uj* . 
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Let us stress that our numerical analysis was performed in Matsubara frequency region, 
which was used in writing down the Eg. (p3D . Analytical continuation of our numerical data 
to the real frequencies was not attempted. 



VI. DENSITY OF STATES CLOSE TO THE METAL-INSULATOR TRANSITION. 

Consider the effects of interelectron interactions upon the single-particle ("tunneling") 
density of states which is defined by the well-known relation: 

N{e) = -lJ-^,ImG^{p,e) (41) 

where e = E — Ep - is electronic energy with respect to the Fermi level and G'-^(p, e) - is 
the retarded Green's function defined by: 

Consider the so called "Fock" contribution to the self-energy part S^(£, p) shown in Fig.6: 

^ z7VG'o (p, e) (/(e, u^) + ig{e, u^)) , (43) 
where /e,^^ and g^^uia defined by the relations: 

A. = 4iV-(0)i?e / ^ ^ , (44) 

g^^ = 4N^\0)Im I ^ C'^ — 2 (45) 

and G(f (p, s) - is the advanced Green's function without interaction contributions. Here and 
in the following we denote as A^o(O) the density of states at the Fermi level in the absence 
of interactions. 

Substituting the expression for S^(£, p) from (|43|) into (^) we obtain as £ ^ 0: 

16 



N{e) , 

iVo(O) TT 



1 °° 

- / dep/mG^(p,£) 

TT J 



-co 



oo 



/ _ _ _ \ 2 



(46) 



^2(1 + f,g,,^^ + ((1 + /i^?,,^J2 + (/i/,,.J2) 1/2) V ((l + /i^?.,.J2 + (/i/,,.J2)l/2^ 

Now let us calculate /^^^^ fl'e,a;c- we noted above, the frequency behavior of diffusion 
coefficient in metallic phase is qualitatively similar to that in the usual self-consistent theory 
of localization (with the shifted mobility edge). In insulating phase, in a narrow frequency 
region u <^ u* the frequency dependence of D^u) is actually unknown to us. However, 
taking into account the fact that as we approach the transition uj* 0, it is reasonable to 
use the expression (5) for the generalized diffusion coefficient D{uj), as given by the self- 
consistent theory of localization. Here we assume that the role of interactions reduces only 
to a simple shift of transition point. 

In metallic region, for e —* we obtain: 



^\EpJ \ uj,'/y n \EfJ XH/ryl^ 



33/2 / 7 \ V \ 33/2 / ^ x2 1/^ 



In the region oi oj^ < e < 1/r the functions f^^uia ^i^d 9e,uic reduce to: 

33 / 7 \3 / ^1/3 



Using the expressions ( ^'^ - |50| ) we can analyze Eq. (|46|) for the density of states N{e). In case 
of relatively weak disorder /^^^^ and gs,uj^ -C 1 and for /i ^ 1 we can restrict ourselves by 
the linear term over fi and obtain the correction found in Refs. [10-11]: 
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"^^^^-■^1-^^...) (51) 



iVo(O) V 2- 

As the system moves towards the mobihty edge (7 ~ Ep, ujc 0) the function (7^^^ diverges 
logarithmically {gs^ujc ~ ^^17") expression (^Tj) is insufficient. We have to 

calculate the density of states using the complete expression (^). Note that in Refs. [10-11] 
the frequency dependence of diffusion coefficient was completely neglected. However, we 
shall see that it will become very important as we approach the transition. For Uc ^ in 
Eg . (^61) we can neglect f^^ujc comparison with the diverging function g^^ui^- Thus we have: 
'VW^,,,,^ 1-1/2 .3"''f7Vf, e^''' . 1 ,„1A~ 



^„(0) - ;i + ' '''-^wAe;)['':F-^ T2"'t) '''' 

In case of the weak disorder the well-known square root singularity appears at the Fermi 
level [3]. As system moves towards the mobility edge {uc — > 0) the divergence of causes 
this minimum to become more deep and at the mobility edge itself {uc = 0) the density of 
states at the Fermi level is equal to zero. Note that the square root dependence of the 
density of states persists only in the region of < e < Uc, and the width of this region tends 
to zero as we move to the mobility edge. In the region of uj^ < e < 1/r we have: 



N{e) / 33/2 / 7 \2 1/ 



-1/2 



while at the mobility edge itself (7 ~ Ep): 



N{e) / 33/2 1/ 



-1/2 



In Fig. 7 (curves 1-3) we show the numerical data for the density of states obtained from 
direct calculations of f^^ujc 9e,uj^ and Eq.(^) which demonstrate the formation and growth 
of effective pseudogap. Curves are given for the different values of parameter A/ Ac, where 
Ac - is the critical disorder parameter determining the metal-insulator transition. Dashed 
curve shows the behavior of the density of states at the mobility edge itself (A/Ac = 1)- 
In insulating phase the functions g^^ujc and f^^uic have the following form as e — 0: 



Neglecting in Eg . (^61) g^^wc comparison with /^^^^ which diverges as e — we get: 

N{e) 1 22 £1/2 



(57) 



iVo(O) /iV2 33/2 ^1/2 

Thus, in the insulating region the square root dependence for e —>■ persists and its region 
widens as disorder grows. For Uc < e < 1/t we obtain the previous dependence given in 
Eq.dD (Fig. 7, curves 4-5). 

Above we have considered the conduction band of infinite width. Important changes for 
the insulating phase appear as we address the case of conduction band of finite width 2. In 
the model of linear electronic spectrum we obtain: 

-(^^-■' f,,,, |±l^ , (58) 



iVo(O) n_J^ (^2+7^) +(/x7V.,.J 
Analysis of expressions (^) and ( p8|) shows that the density of states in the region of 
£ <^ fxcUc demonstrates quadratic e - behavior: 

N{e) f bV [EfV 1 / e 



NoiO) \lj V 7 ; f^' ' ^^^^ 

which reminds the well known "Coulomb gap" of Efros and Shklovskii [13-14] valid deep 
within insulating phase. Note that this may a pure coincidence because the analysis of Efros 
and Shklovskii is essentially based upon the long-range nature of Coulomb interaction. The 
numerical results for N{e) are shown in Fig. 8 (curves 5-6). In metallic region in the case 
of conduction band of finite width numerical data demonstrate that the density of states is 
practically unchanged in comparison with the case of the infinite band. 

Tunneling densities of states dependences similar to those shown in Fig. 8 were observed 
experimentally for a number of disordered systems close to the metal-insulator transition 
[1,2], from amorphous alloys [15-17] to disordered single- crystals of metallic oxides, including 
high-temperature superconductors [18]. While we see the complete qualitative correspon- 
dence of theory and experiments, though we must stress that the quantitative agreement is 
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obviously absent — the effective width of the "pseudogap" close to the transition is larger 
in the experiments and can be described by power-like, not by logarithmic dependence. At 
the same time our expressions describe the significant overall drop of the density of states in 
rather wide energy region, which is observed in all experiments. The expressions obtained 
above give, as far as we know, the first description of the evolution of the tunneling density 
of states during the transition from metallic to insulating state induced by disordering, for 
the whole interval of parameters controlling this transition. 

VII. CONCLUSIONS. 

In this paper we propose a self- consistency scheme which allows the description of disorder 
induced metal-insulator transition taking into account the first-order effects of interelectron 
interactions. This approach is based upon the assumption that neglecting the interelectron 
interactions this transition can be sufficiently well described by self-consistent theory of lo- 
calization, while interaction effects may be taken into account if we consider the appropriate 
interaction and localization diagrams on the equal footing in the equation for the general- 
ized diffusion coefficient. We demonstrate that this approach allows us to formulate unified 
equations describing the transition from metallic to insulating state and obtain a number 
of results which are in qualitative agreement with experiments on disordered systems (tun- 
neling density of states). Interaction effects remain relatively small in this approach and 
for the wide region of parameters the diffusion coefficient is well described by self-consistent 
theory of localization. In particular, metallic conductivity (for T — 0) linearly drops to zero 
at the transition, which is observed in a number of real systems [1,2] However, this approach 
is unable to describe the known class of disordered systems where this conductivity drop is 
apparently described by the square root dependence on disorder parameter [1,2]. 

From purely theoretical point of view, even if we restrict ourselves only to first-order 
interaction corrections, we still need to perform a complete analysis of "Hartree"-type in- 
teraction corrections neglected above, as well as the importance of screening anomalies in 
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insulating phase [5,11]. We have also neglected all the effects due to electronic spin, while 
the growing importance of these effects as we approach the metal-insulator transition were 
noted long ago [1,2]. 

Still, the proposed scheme of analysis is, in our opinion, of considerable interest because 
it allows us to take into account interaction effects in a wide region of parameters controlling 
the metal-insulator transition, including even the insulating phase, which was inaccessible 
for the previous theories. 

This work was partially supported by Russian Foundation for Fundamental Research 
under the grant N° 93-02-2066 as well as by Soros International Science Foundation grant 
RGLOOO. 
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Figure Captions. 



Fig.l. Two-particle Green's function $pp/(qc<j£:) and vertex parts 7^^(pqc<j) and 
Fig. 2. First-order interaction corrections to conductivity. 

Fig. 3. Static generalized diffusion coefficient (d = -75^) dependence on disorder in dif- 
ferent self-consistency schemes: (a) — scheme I; (b) — II; (c) — III; (d) — IV. 

Fig. 4. Dependence of dimensionless generalized diffusion coefficient on dimensionless 
Matsubara frequency in metallic phase (a = 0.5), obtained by numerical solution of Eq.(p3D 
for different values of fi: 1. 0,24; 2. 0,6; 3. 0,95; Dashed line — the usual self-consistent 
theory of localization, fi = 0. 

At the insert: Dependence of static diffusion coefficient {d = on disorder for fi = 
0,24. 

Fig. 5. Dependence of dimensionless generalized diffusion coefficient on dimensionless 
Matsubara frequency in dielectric phase {a = —0.5), obtained by numerical solution of 
Eq.(p3D for different values of fi: 1. 0,12; 2. 0,6; 3. 1,2; Dashed line — the usual self- 
consistent theory of localization, fi = 0. 

Fig. 6. "Fock" contribution to self-energy part. 

Fig. 7. Evolution of the density of states during the metal-insulator transition in the 
interacting system with /i = 0.2. Infinitely wide conduction band. Data are shown for 
different values of disorder parameter A/Ac: 1.0.5; 2.0.7; 3.1; 4.1.4; 5.1.8 

Fig. 8. Evolution of the density of states during metal-insulator transition in the inter- 
acting system with fi = 0.3. Conduction band of finite width 2B = AEp. Data are shown 
for different values of disorder parameter A/Ac: 1.0.5; 2.0.7; 3.0.9; 4.1; 5.2; 6.2.4 
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